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Dynkin’s {Bull. Amer. Math. Soc. 3 (1980) 975-999) seminal work associates a multidimensional 
transient symmetric Markov process with a multidimensional Gaussian random field. This as¬ 
sociation, known as Dynkin’s isomorphism, has profoundly influenced the studies of Markov 
properties of generalized Gaussian random fields. Extending Dykin’s isomorphism, we study 
here a particular generalized Gaussian Markov random field, namely, the de Wijs process that 
originated in Georges Matheron’s pioneering work on mining geostatistics and, following McGul- 
lagh {Ann. Statist. 30 (2002) 1225-1310), is now receiving renewed attention in spatial statistics. 
This extension of Dynkin’s theory associates the de Wijs process with the (recurrent) Brownian 
motion on the two dimensional plane, grants us further insight into Matheron’s kriging formula 
for the de Wijs process and highlights previously unexplored relationships of the central Markov 
models in spatial statistics with Markov processes on the plane. 
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1. Introduction 

After originating in the pioneering work of Georges Matheron, the de Wijs process enjoyed 
a significant and extensive role in early geostatistical literature [6, 19, 26, 28]. McCullagh’s 
[30] recent work has revived interest in the de Wijs process, both theoretically and in a 
growing range of applications in spatial statistics; see, for example, [4, 5, 8-10, 13, 31, 34- 
36]. In particular, Mondal [34] and Besag and Mondal [5] established a connection between 
Gaussian Markov random fields on two-dimensional lattices and the de Wijs process on 
the Euclidean plane, which emerges as a scaling limit of the former. McCullagh and 
Clifford [31] analyzed agricultural uniformity trials using a spatial formulation that is 
based on the de Wijs process and a Gaussian white noise random held. See also the related 
work by Clifford [8, 9], Clifford et al. [10]. Mondal [35] considers the exponential functional 
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of the de Wijs process to construct a generalized Cox process to study disease mappings. 
Mondal [36] indicates a link between the de Wijs process and Tobler’s [52] pycnoplylectic 
interpolation based on the Laplace splines. Dutta and Mondal [13] make explicit use of 
the connection between intrinsic autoregressions and the de Wijs process and provide 
approximate matrix free computations for residual maximum likelihood methods for the 
latter. Furthermore, outside the statistics literature, the de Wijs process appears to 
originate separately in quantum physics and statistical mechanics as the massless case 
of the free Gaussian field; see, for example. Chapters 6 and 7 of [17], and in recent 
probability literature, this massless case has become a subject of intense study; see, for 
example, [21, 48]. 

Technically, the de Wijs process is a generalized Gaussian random field ([16], Chapter 
111), whose index set is a certain class of contrasts, that is, non-atomic signed Borel 
measures on the Euclidean plane with zero total mass. This process corresponds to the 
logarithmic variogram model and is a generalization of the Brownian motion in two 
dimensions. It acquires Markov and conformal invariance properties [30, 31] and is first- 
order intrinsic in the sense of [55] and [12]. The Markov property of the de Wijs process 
was already known to Matheron [27, 28], who viewed it from the perspective of kriging 
predictions. Consider a mean zero Gaussian random field {U{x) :x € TZ^}. Let B denote 
a closed contour in TZ?. It is natural to call the random field Markov if its values along 
the curve B determine the ordinary kriging predictor for the value U{xq) at a point 
xq in the interior of B, given its values on and in the exterior of B. This Markovian 
characterization leads to a kriging predictor for the random variable U{xo) that takes 
the form of a contour integral 



where the coefficient function w(a;,a;o), x S i? is such that 



( 1 . 1 ) 


for every point x' on B. When {U{x ): x G TZ?} is stationary and isotropic, Matheron [27] 
deduces that it is Markov if and only if 


cov{U{x),U{0)) ociLo(alkll), 


( 1 . 2 ) 


where Kq is the Bessel function of order zero and a is a positive constant. Matheron 
further notes that his derivations remain intact for the limiting case a j, 0 that corre¬ 
sponds to the logarithmic covariance, cov(f/(x),t/(0)) = — log(||xll) and thus provides 
the Markovian characterization of the de Wijs process. However, the respective random 
fields exist as generalized processes only, and the details of a formal argument go beyond 
the above kriging formula. 

This paper calls attention to the work of Dynkin [14] to present a mathematical for¬ 
malism to describe the above kriging formula of the de Wijs process, and to connect the 
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field of spatial statistics to the vast, and hitherto unutilized, probabilistic literature on 
the Markov property of generalized random fields. This body of literature constitutes a 
fascinating part of probability, much of which emerged in the wake of [24] and [32] . This 
corpus notably includes [i4, 15, 18, 20, 22, 25, 33, 38, 39, 42-45, 47, 53, 54] and the ref¬ 
erences therein. In many of these works, several notions of Markovianity for generalized 
Gaussian random fields have emerged and their interrelations and their connections to 
various related concepts often form a good part of their understanding. For example, a 
homogeneous and isotropic generalized Gaussian random field whose spectral density is 
inversely proportional to an even polynomial of the frequencies satisfies a Markov prop¬ 
erty in the sense of Holley and Stroock but may not be Markov in the sense of Wong. It 
is interesting to note that Nelson’s [38] construction of the free Markov field on the plane 
actually corresponds to the generalized Gaussian random field with covariance given in 
(1.2). Using a slightly different notion of Markovianity, Wong [54] arrives, much earlier 
than Nelson, at the conclusion that the only generalized Gaussian Markov random held 
again has covariance (1.2). Kallianpur and Mandrekar [20], on the other hand, investigate 
Markov properties of a generalized Gaussian random held in conjunction with its dual 
random held. Ekhaguere [15] later provides links between the Markov property due to 
Nelson [38] and that due to Wong [54]. In contrast, Dynkin’s [14] famous work marks an 
important departure from these earlier studies. In his study, covariances of a generalized 
Gaussian random held are assumed to arise from the Green function of a symmetric 
multidimensional Markov process, and the Markov property of this generalized Gaussian 
random held is then derived from the path properties of the multidimensional Markov 
process. Thus, for example, the Markov property of (1.2) can be understood in the con¬ 
text of the Markov property of an exponentially killed Brownian motion on the plane. 
Here our focus is on the limiting case, namely, the de Wijs process. Although its Markov 
property can be investigated using the work of Nelson or Wong (e.g., by modifying The¬ 
orem 1.5 of [18] or by including the case a = 0 in Wong’s [54] Theorem 2), we take up 
Dynkin’s approach primarily because it provides a precise and computable description 
of the boundary condition in the kriging formula, and connects closely with Matheron’s 
work. We also piece together many scattered results and extend some known ones to 
provide this new addition to the body of literature that, respectively, followed Dynkin’s 
and Matheron’s works. 

The remainder of the paper is structured as follows. Section 2 introduces the de Wijs 
process as a homogeneous, isotropic and self-similar generalized Gaussian random held. 
Section 3 explores the association of this process with Brownian motion. Here we show 
that the covariance formula of the de Wijs process can be written explicitly in terms of 
an additive function of the Brownian motion. Section 4 studies the Markov property of 
the de Wijs process by extending the work of Dynkin [14]. Here our main result, namely 
Theorem 4.1, provides a new interpretation of Matheron’s kriging formula in terms of 
the hitting probabilities of the Brownian motion and as a generalization of the Dirichlet 
problem. Section 5 focuses on the practical relevance of Matheron’s kriging formula. It 
also considers the relevance of Dynkyn’s isomorphisms in lattice approximations of the 
de Wijs process and concludes with a discussion on the screening effect in kriging. 
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2. De Wijs process 

In this paper, a generalized random field on the Euclidean plane Ti? is a stochastic process 
{Zcr-aGM} indexed by a vector space A4 of non-atomic signed Borel measures on the 
plane that have total mass zero. We view Z to be a linear functional from the vector 
space Ai to the real numbers such that 

Zba+di' = bZc + dZjy for all a,v ^ AA., and for all b,dGTZ. 

We think of the random variable Z^ as a spatial contrast \ for instance, if two plots have 
unit area and tr has a Lebesgue density that is proportional to the difference of the 
respective indicator functions, then Z„ might represent the difference of crop yields on 
these plots. 

The generalized random field {Z^:aGM.} is said to be homogeneous if its distribution 
remains invariant to planar translations, and isotropie if its distribution remains invariant 
to planar rotations. Furthermore, such a homogeneous isotropic generalized random field 
is Gaussian if all finite dimensional marginal distributions are multivariate normal with 
EZo- = 0 and 

cov{Z^,Z^) = - JJ ip{\\x - y\\)a{dx)v{dy) (2.1) 

for some real-valued function ip and all non-atomic signed measures a,v € AA. Note that 
in the above covariance formula we can add to (/9(||x —?/||) a function fi{x) + f 2 {y) without 
affecting the integral, and so ip actually belongs to a suitable quotient space (modulo the 
infinite dimensional subspace of additive functions). In subsequent discussions, we will 
be implicit about this equivalence relation in the description of ip. The function 

c(x) = 1 ^( 11 x 11 ), 

is then called the generalized variogram or —c{x) the generalized covariance function of 
the generalized random field. Let d and z) denote the Fourier transforms of a and v. We 
can then write (2.1) as 

cov{Z„,Zi,)= j d{x)v{x)S{dx) (2.2) 

for a certain non-negative tempered measure S which is called the spectral measure 
of the generalized random field ([16], page 264). If the spectral measure is absolutely 
continuous with respect to the Lebesgue measure on TZ^, its Lebesgue density s(x), 
X G TZA is called the spectral density. Under slight regularity conditions, the generalized 
covariance function c and the spectral density s are Fourier transforms of each other. 

Specifically, consider Ad to be the space of signed Borel measures a on the Euclidean 
plane VA that satisfy 


|log(||a;- 2 /||)||cr|(dx)|cr|(d 2 /) < oo 


(2.3) 
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and have total mass zero. The de Wijs process is then the homogeneous, isotropic and 
self-similar generalized Gaussian random field (i.e., its distribution also remains invariant 
to changes of scale) on TZ^ with index set M such that EZ^ = 0 and 

cov(Z^, = {a, v)m = - yy log(||x - y\\)a{dx)v{dy) (2.4) 

for all signed measures a,v &M. Note that M has an inner product space structure with 
inner product (cr, v)m a-nd norm 


\W\\m 


(cr, cr) 


1/2 

M ■ 


Indeed, Corollary 2.5 of [29] implies that M is a. vector space, and by Corollary 2.4 
and Remark 3.3 in the same reference ||cr||^ >0 with equality if and only if cr vanishes 
identically. The positive definiteness of the covariance matrices associated with the de 
Wijs structure (2.4) is an immediate consequence of the Gram matrix property. Thus, the 
de Wijs process has logarithmic variogram; that is, the representation (2.1) holds with 
the generalized variogram function c{x) = (/5(||x||) = log ||a;|| for x G 7^^, and its spectral 
density is 


s(a;) 


1 

2 ^’ 


xGTZ^. 


3. Association with Brownian motion 

We set T = [0,oo) for consistency in what follows and let {Wt,t G T} be the Brownian 
motion on the two-dimensional Euclidean plane. Thus, with probability 1, the function 
t ^W{t) is continuous in t, the components of the increment Wt+u — Wu are independent 
Gaussian random variables each with mean 0 and variance t, and the process {Wt,t G T} 
has stationary and independent increments. For every x on the plane, let P^, denote the 
probability law of {Wt,t G T} starting at x and let E^, be its expectation under P^,. 
For every t gT, let the sub-tr-field consist of events observable up to time t, which 
is the minimum cr-field generated by {Wu :0 < u < t}. Define to be the minimum 
O'- field containing The Markov property of the Brownian motion asserts that 

the conditional law of {Wt,t > u} given {Wt, 0<t<u} depends on W only through Wu- 
In other words, 

E,(FJ) = E,(EEw„J) 

for every x on the plane, every ^u measurable positive function F and every measurable 
function J that depend only on {Wt^t > u}. In particular, the expectation is calculated 
first with respect to the conditional law of {Wt,t > u} given {Wt,0 <t< u}, and then 
with respect to the marginal law of {Wt,0 <t<u}. An important generalization of 
the Markov property is the strong Markov property. When r is a stopping time, define 
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the stopping field to be the cr-algebra consisting of all events A G such that 
^ n {r < t} G for every t > 0. Then the strong Markov property implies that 


E,(FJ)=E,(FEiv^J) 

for every x on the plane, every stopping time r, every measurable positive function 
F and every measurable function J that depends on {WtA > r}. 

Next we define the Green function of the Brownian motion {Wt,t € T}. Typically, 
the Green function is defined for a transient Markov process as the time integral of its 
transition probability density function, and Dynkin’s theory is essentially based on the 
fact that the Green function of a transient symmetric Markov process can be interpreted 
as the covariance of a centered Gaussian process. However, the Brownian motion on the 
plane is recurrent and hence its transition probability density function 


Vt{x,y) = {2nty 


exp<{ -^||y-x|| 


is not integrable with respect to t. Thus, we need a modification that will allow us 
to define the Green function of the Brownian motion and extend Dynkin’s result in a 
straightforward fashion. To this end, we fix a point xq on the unit circle and consider 
qt{x,y) = pt{0,y — x) —pt{0,xo)- We then apply the definition of Port and Stone ([40], 
page 70) and obtain the Green function or the potential kernel of {Wt,t G T} as 


poo 

gix,y)= / qt{x,y) dt = -7X.-'^ log\\y-x\\. 
Jo 


(3.1) 


The choice of xq will not matter, as we shall see in what follows. 

However, note that the covariances of the de Wijs process now satisfy the relationship 

cov{Zc,,Zy = {a,v)M = j j T^g[x,y)a{dx)v{dy). 

Thus, we say that the above relationship associates the de Wijs process with the Brownian 
motion {Wt,t G T}, opening an avenue for exploring the properties of the former from 
those of the latter. For every z/ G A4, we now get 

J qt{x,y)iyidy) = J pt{x,y)iy{dy) - pt{0,\\xo\\) J v{dy) = J pt{x,y)v{dy) 


and therefore 


j 9{x,y)i^{dy)= Jj pt{x,y)v{dy)dt. 


In other words, the term involving xq disappears from the right-hand side of the previ¬ 
ous equation. Consequently, when v is absolutely continuous with the Radon-Nykodyn 
derivative v(dy) = p{y) dy, the above equation becomes 


g{x,y)iy{dy) = E^ J p(Wt)dt. 
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Now equation (3.1) can be identified with J p{Wt)dt, where Ecr is the 

expectation under the probability law of {Wt,t gT} with initial signed measure tr; that 
is, 

Ect J piWt)dt = J J p(Wt)dt a(dx). 

We can thus define the additive function of the Brownian motion by the measure 

A,iQ)= [ p(Wt)dt (3.2) 

Jq 

that satisfies the property that, for every interval I = {s,u) with s <u, A^(I) is a func¬ 
tional of G /}, and 


cov{Z^,Z^) = =nEaA^{T). (3.3) 

The collection of all signed measures v € A4 that are absolutely continuous forms a dense 
subspace of A4. By passage to limit, it then follows that for every i/ € Ai there exists 
an additive functional of the Brownian motion such that the above equation holds. In 
addition, the strong Markov property of the Brownian motion takes the following form 

E„FA,{t + Q)=E,FEwMQ) ( 3 - 4 ) 

for every cr, i/ G A4, for every Borel subset Q of T, and for every r and F as defined 
earlier. The strengthened relationship that emerges from equation (3.3) in conjunction 
with equation (3.4) now paves the way to use the Brownian paths to study the properties 
of the de Wijs process. 


4. Markov property of the de Wijs process 


As in Section 1, take S to be a simple closed contour on the plane. Then B divides 
the entire plane into two components, namely, the bounded interior and the unbounded 
exterior. Let Bi denote the open interior of B with closure Bi. Similarly, let Be be the 
open exterior of B with closure Be- Our first task is to describe the values of the de Wijs 
process on the boundary B, and on the inside and the outside of B (e.g., on sets Bj and 
Be)- To this end, there are two approaches. The first approach is due to [39]. Here we 
describe the values of the de Wijs process on an open set G by the minimum sigma field 
Ag generated by all Z^ such that cr G A4 and support of a is compactly contained in G. 
Then, for any closed set G the values of the de Wijs process is described by the sigma 
field 


Ac = Pi Ag, 

gdc 

where the intersection is taken over all open sets G that contain G. Thus, Abi, -^Bb 
and Ab represent the values of the de Wijs process on the inside, outside and on the 
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boundary, respectively, and the Markov property of the de Wijs process asserts that, for 
any a € M with support of a compactly contained in Bi, 


E{ZMsJ=E{Z,\Ab) 

almost surely in the probability distribution of the de Wijs process. Note that, by con¬ 
struction, the minimum sigma fields and Ab contains neighborhood information, 
not just the information on the set. The second approach adopted by Dynkin [14] is a 
simplified version of the above and goes as follows. For a close set C, we define A4c to 
be the set of signed Borel measures a € A4 that do not charge on its complement. Then, 
following [14] , the minimum sigma field generated by the collection of random variables 
{Za : cr G A4c} describes the values of the de Wijs process on the C. The Markov prop¬ 
erty is described in the usual fashion, namely, the values on the inside (Bi) and the 
outside (Be) of B are conditionally independent given the values on B. This leads to the 
following theorem. 


Theorem 4.1. Let D he any closed set on the plane with a simply connected open 
interior, and let tb he the first hitting time of D by the Brownian motion {Wt,t G T} 
starting at x on the plane. Denote by 14 the probability measure of IFVj, conditioned 
on Wq = X. For any initial signed measure a € A4 which does not charge on the set D, 
identify ud with the signed measure induced by Wt at the first hitting time of D; that is, 

aD{G)= Jp,(Wrr,GG)a(dx) = J V,,(G)a{dx) 

for all Borel subsets G on the plane. Then the conditional expectation of Z„, given the 
values of the de Wijs process on D, is identical to . In other words, 

E(Z,\{Z,,v^MDf) = Z,^. (4.1) 


Before we turn to the proof, let us first see how Matheron’s kriging formula enters into 
the above theorem. Equation (4.1) implies that for all tr G A4 and v & M.d 


EZ^Zi, — EZo-jjZy. 


Hence, formula (2.4) applies. This along with the definition of ao produces the identity 

J j logdlx - y\\)a(dx)iy{dy) = JJ J logdk' - y\\)V^(dx')(7{dx)v(dy) 

for all tr G A4 and v G JAd- Consequently, if x is in the interior of boundary B, we have 

[ V:^(dx') = l, log(||x- 2 /||)= [ logdla;'-j/ll)14(da;') 

J B J B 


for every point y on B. In short, the coefficient function v(x',xo) in Matheron’s kriging 
formula (1.1) is the derivative of t4o A x', and thus corresponds to the probability 
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density function of the Brownian motion at the first hitting time td, a crucial fact that 
has arguably been missing from the geostatistitical literature. Furthermore, in order for 
Z to be a linear functional from the vector space A4 to the real numbers, we can imagine 
as an integral of the form J ZxCr(dx), where the notation Z^ suggests a point-wise 
intrinsic process with vai(Zx — Z^') = — logdlx — x'H). This very imagination of a point- 
wise Zx allows us to describe Z^ = /g Zx'Vx(dx') as the kriged value of Z^, for an x in 
the interior of B. We can take this point further, and even describe the kriging formula 
from a different angle. First, let V denote the Laplace operator on the plane. If u is twice 
differentiable. Theorem 3 of [41], page 525, implies 

— 2Ti{(7v)M = J cF{x)iy{x) dx, (4.2) 

which surprisingly asserts that Z^ and Z^' behave as two mean zero uncorrelated 
Gaussian random variables for x ^ x', and in turn, suggests that the kriged values of 
VZj, on the interior of B are all zero. An interchange of the Laplace operator and the 
conditional expectation on Z^ then produce 

Zx = 0 

on the interior of B, implying that the kriging problem is a generalization of the Dirichlet 
problem in mathematics. Indeed, the probability literature reaffirms that the boundary 
values of the Brownian motion at first hitting time solve the standard version of the 
Dirichlet problem and therefore the Matheron’s kriging formula (1.1) can be seen as a 
generalization. For significance of the Dirichlet problem in recent spatial statistics, we 
refer to the discussion in [30]. We now return to the proof of the theorem. 

Proof of Theorem 4.1. First, we verify that ao € Md- As a first step, we argue 
that an belongs to M. Since f a^idx) = f (T(da;) = 0, an represents a signed Borel 
measure with total mass zero. Now, for a non-negative measure p, for which the integral 
hfj,{x) = / g{x,y)fi{dy) is finite for every x, the results of [7], pages 193-194, give the 
identity 

Ps{x,y)ds p.{dy) > 0. 

Consequently, defines an excessive measure, and h^(x) > Bix:h^{fWt)- Hence, the choices 
pL = v^ and g = v~ yield 

K+{x) > Exh,,+ {Wt), K-{x) > ¥,x:K-{Wt). 

It then follows that 


hfj,{x)-E,xhfj,{Wt)= / g{x,y)p.{dy) - Ex / giWt,y)g{dy) = 


h,,+ {x)(T+{dx) = {a+,v+)j^, 
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and, after repeating the same argument, )m < (c ,v )m and so on. Thus, 

^ j^\ + \{(j+ ,v-) 

and a similar upper bound exists for {<j^,v)m- The above bounds imply 

{(JD,(Jd)m < + 2|(cr+,cr")yKl + I (o’” , cr”)^ |, 

which ensures that an G A4. Now H is a closed set with a simply connected open interior 
and so we get 

|ar3|(L»“) = J ^i?)|a|(dx) = 0. 

Therefore, an belongs to A4d- Next, we establish that 

= ^i'&A4n- 

Since an is a measure that satishes the relation an{G) = Y^ai^dWro))^ for all Borel sub¬ 
sets G of the plane, integrals with respect to an can be defined as appropriate expected 
values of the functions of Wr^ ■ In particular, for any element / of an appropriate class 
of functions, such an integral will satisfy 

J f{x)an{dx)=F.^f{Wro)- 

Now take f{x) = FixA,^{T). Then, the definition of the additive function asserts that 
FjZo-jjZi, = nEcr^A^(T), but the above equation also implies 

E,^A,(T) = J EMT)anidx)= J f{x)anidx)=E,fiWr^)=E,Ew^^A,iT). 

Consequently, the strong Markov property of the Brownian motion in equation (3.4) 
applies, and we obtain 

EacAuiT) = E^Ayijn + T). 

Since Tn is the first hitting time of D, the path of the Brownian motion up to but not 
including time Tn lies entirely within the complement of D. However, the signed measure 
V concentrates on D making it imminent that E^Ai,((0,Tn)) = 0. And, therefore 

cov{Zcro,Z^) =nEcroA^{T) =7zE^A^{Tn +T) =7iEcr A^{T) =cov{Z^,Z^). 

This completes the proof. □ 

Interestingly, when values are known along a straight line or on a circle, the ana¬ 
lytic formulas for the coefficient function v(x,xo) are available in closed form, making it 
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possible to apply Theorem 4.1 directly to calculate relevant kriging predictions. For an 
example, when B is the unit circle, v{x,X{f) becomes the Poisson kernel 

Furthermore, when the boundary set B is the y-axis, we refer to [27] for a formula for 
the corresponding coefficient function v{x,xq). 

Finally, we can also discuss predictions for functionals of the values of the de Wijs 
process inside the boundary B (e.g., f{Z„) for some suitable function /), given {Z^ : v £ 
Md}, but this would require a knowledge of Wick products and Fock spaces and is 
beyond the scope of this paper. 

5. Discussion 

In practice, we only select finitely many regular or irregularly distributed sampling lo¬ 
cations and observe process values as aggregates or averages over certain non-empty 
regular or irregular regions around those sampling locations. In both instances there are 
certain limitations in applying Matheron’s kriging formula directly. For example, if we 
observe only finitely many data values on the unit circle, we won’t be able to apply the 
exact kriging formula in (4.3). Similarly, when the de Wijs process is used as a statisti¬ 
cal model for aggregates or averages of spatial variables over non-empty regions in the 
two-dimensional plane, we simply lose the Markov property because of aggregations or 
averaging. Examples include agricultural field trials where the variable of interest is the 
crop yield over plots, or disease mapping where the spatial variable of interest is consid¬ 
ered a stochastically degraded version of an underlying unobserved spatial component 
such as the log relative risk of non-infectious diseases over a geographic region. 

However, certain discrete approximations are possible, and, in fact, it is the discrete 
approximations of the de Wijs process that have played a major role in spatial statistics in 
the past thirty years; see, for example, [2, 3, 23], many subsequent papers, and the books 
by Cressie [II], Banerjee et al. [I] and Rue and Held [46]. These discrete approximations 
form a subclass of Gaussian Markov random fields on regular and irregular lattices and 
have lattice graph Laplacians as their precision (i.e., inverse of ‘covariance’) matrix. The 
nature of these approximations become clearer when we also note that the log function 
(i.e., the generalized covariance of the de Wijs process) is the inverse of the Laplacian 
on the plane. Interestingly, these Gaussian random fields are also associated with the 
random walks on the lattice graph, as the De Wijs process is with the Browning motion 
on the plane. As a concrete example, the first order symmetric intrinsic autoregressions 
on the two dimensional integer lattice [5, 34] is associated with the simple random 
walk on Z^. Thus Dynkin’s theory also applies here and we can obtain the coefficient 
function of a corresponding kriging problem on the discrete lattice Z^ from probabili¬ 
ties of the simple random walk at the first hitting time. To summarize, the diagram in 
Figure 1 lists a few important spatial Gaussian Markov models and the associated two- 
dimensional Markov processes. The top part of this diagram notes the lattice Gaussian 
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Stationary autoregression on Intrinsic autoregression Z^ 

s{uj,7j) = [1 -^ + /3{sin^(ia;) + sin^(i77)}]"^ ..> s{uj,7j) = (sin^(ia;) + sin^(ir?))“^ 

(Simple random walk with (Simple random walk) 

geometric holding times) 

1 

Geireralized Ornstein-Uhlenbeck process 
s{iJ,r]) = (a^ + ??^)“^ 

(Brownian motion with 
exponential holding times) 


De Wijs process 
s{uj,r,) = {u^ + r,^)-^ 
(Brownian motion) 


Figure 1. Limit diagram for Gaussian Markov random fields and associated Markov processes. 
Here 0 < /3 < 1 and a > 0. 


Markov random fields along with the spectral densities and, in brackets, the associated 
lattice Markov processes. The bottom part of the diagram provides limiting continuum 
Gaussian Markov random fields with corresponding spectral densities and, in brackets, 
the associated Markov processes. These continuum random fields arise as the scaling 
limits of corresponding lattice Markov fields from the top part of the diagram; see [37] 
and [5] for details. 

Another interesting point is that aggregates or averages of the de Wijs process retain an 
approximate Markov property that is known as the screening effect in geostatistics [6, 50]. 
To give a very simple example, we consider a regular lattice in the two-dimensional plane. 
Let Xg^t denote the average value of the spatial variable of interest over the unit square 
whose center has integer coordinate (s, t) in the Euclidean plane. Given a realization 
of Xg^t for s,t = —8, ...,8 but excluding Xo,o, we can employ exact variogram 
computations [8, 34] to find the coefficients ujg^t in the conditional expectation or ordinary 
kriging predictor [49] 

( 5 - 1 ) 

under the regularized de Wijs process. The sum on the right-hand extends over the 
aforementioned index set and the ordinary kriging coefficients ujg^t add up to 1. Table 1 
shows the numerical values of cog^t to three decimals; in view of symmetries only 44 of 
the 17^ — 1 = 288 coefficients need to be shown. The screening effect is prominent here 
in that the immediately neighboring cells dominate, with very few of the remaining cells 
receiving non-negligible ordinary kriging coefficients. However, it is not known to me if 
uig^t can be interpreted in a meaningful way using certain probability calculations of the 
Brownian motions, but here one can further try to derive analytic form of cog^t for an 
infinite lattice from the spectral density form of X^^v The same applies for aggregates 
or averages of the first-order intrinsic autoregression. In general, it would be interesting 
to know if one can better understand such an approximate Markov property. 
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Table 1. Numerical values of the coefficients ojs.t in the ordinary kriging predictor (5.1) under 
the regularized de Wijs process 


s 

t 








1 

2 

3 

4 

5 

6 

7 

8 

0 

0.342 

-0.075 

0.017 

-0.004 

0.001 

0.000 

0.000 

0.000 

1 

-0.032 

-0.001 

0.002 

-0.001 

0.000 

0.000 

0.000 

0.000 

2 


0.002 

-0.001 

0.000 

0.000 

0.000 

0.000 

0.000 

3 



0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

4 




0.000 

0.000 

0.000 

0.000 

0.000 

5 





0.000 

0.000 

0.000 

0.000 

6 






0.000 

0.000 

0.000 

7 







0.000 

0.000 

8 








0.000 


Some future directions can be added to this work. For example, the work of [51] 
provides links between Dynkin’s isomorphisms and constructions of statistical designs. 
Generalizations of Tjur’s work in the context of spatial designs would be an interesting 
matter for future study. 
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